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Considerations of local atomic level stresses associated with each atom represent a particular 
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correlations in a two-dimensional model liquid using molecular dynamics simulations in the following 
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[Phys. Rev. E 91 , 032301 (2015)] related not to the alleged anisotropies of the Eshelby’s stress 
fields, but to the rotational properties of the stress tensors. 
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I. INTRODUCTION 


It is relatively easy to describe structures of crystalline 
materials due to the presence of translational periodicity. 
This periodicity implies that atoms whose coordinates 
differ by a vector of translation have identical atomic 
environments. In glasses and liquids, in contrast, ev¬ 
ery atom, in principle, has a unique atomic environment 
[1, 2]. Largely for this reason description of disordered 
materials continues to be a challenge. Many different 
approaches have been suggested to describe disordered 
structures. However, none of them allows establishing a 
clear link between the structural and dynamic properties 
of disordered matter [1, 2]. 

The concept of local atomic level stresses was intro¬ 
duced to describe model structures of metallic glasses 
and their liquids [2-6]. For a particle i surrounded by 
particles j, with which it interacts through pair poten¬ 
tial U{rij), the aP component of the atomic level stress 
tensor on atom i is defined as [3-6]: 
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The sum over j in (1) is over all particles with which 
particle i interacts. In (1) V) is the local atomic volume. 
By convention, the definition without V) corresponds to 
the local atomic level stress element [7, 8]. Note that a- 
component of the force acting on particle i from particle j 
is /“ = [dU{rij)/dr,j] (rfjirij), where = fj-n is the 
radius vector from i to j. Also note that the atomic level 
stress tensor (I) is symmetric with respect to the indexes 
a and /3. Thus in 3D it has 6 independent components 
[4-6], while in 2D it has 3 independent components. 


There are several important results associated with 
the concept of atomic level stresses. One result is the 
equipartition of the atomic level stress energies in liquids 
[4-6]. Thus the energies of the atomic level stress com¬ 
ponents were defined and it was demonstrated for the 
studied model liquid systems in 3D that the energy of 
every stress component is equal to kifT/A. Thus the total 
stress energy, which is the sum of the energies of all six 
components, is equal to 6-kbT/A = {3/2)kbT, i.e., the po¬ 
tential energy of a classical 3D harmonic oscillator. An 
explanation for this result has been suggested [4-6]. The 
equipartition breaks down in the glass state. Then there 
was an attempt to describe glass transition and fragili¬ 
ties of liquids on the basis of atomic level stresses [9]. 
Another result is related to the Green-Kubo expression 
for viscosity. Thus the correlation function between the 
macroscopic stresses that enter into the Green-Kubo ex¬ 
pression for viscosity was decomposed into the correlation 
functions between the atomic level stress elements. Gon- 
siderations of the obtained atomic level correlation func¬ 
tions allowed demonstration of the relation between the 
propagation and dissipation of shear waves and viscos¬ 
ity. This result, after all, is not surprising in view of the 
existing generalized hydrodynamics and mode-coupling 
theories [10, 11]. However, in Ref. [7, 8, 12, 13] the issue 
has been addressed from a new perspective and the rela¬ 
tion between viscosity and shear waves was demonstrated 
very explicitly. 

Recently it has been claimed in Ref. [14] that considera¬ 
tions of the correlations between the atomic level stresses 
allow observation of the angular dependent stress fields 
which are present in liquids in the absence of any exter¬ 
nal shear. In many respects our attempt to understand 
the results presented in Ref. [14] lead to the present pub- 
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lication. 

Here we demonstrate that the angular dependencies 
of the stress correlation functions presented in Ref. [14] 
do not correspond to the angular dependent stress fields 
(which can exist in the system). We show that the an¬ 
gular dependencies of the stress correlation functions ob¬ 
served in Ref. [14] originate from the rotational properties 
of the stress tensors. 

However, the ideas presented here go beyond the scope 
of Ref. [14]. Here we address the atomic level stresses and 
correlations between the atomic level stresses of atoms 
separated by some distance from a new and yet very natu¬ 
ral perspective. It is surprising that this approach has not 
been investigated in detail before. Reasoning in a similar 
direction was presented in Ref.[15, 16]. However, con¬ 
siderations presented there do not address correlations 
between the atomic level stresses of different atoms. 

This paper is organized as follows. In section H the 
idea of the approach is presented. Section HI is a re¬ 
minder about transformational properties of the stress 
tensors. Atomic level stress correlations functions are dis¬ 
cussed in the context of the present approach in section 
IV. In section V the connection between the Eshelby’s 
inclusion problem and the atomic level stress correlation 
functions is analysed. In section VI the results of our MD 
simulations are described. We conclude in section VH. 


II. STRESS TENSOR ELLIPSES 

The atomic level stress tensor defined with equa¬ 
tion (I) is real and symmetric. Thus it can be diago- 
nalised and, in 2D, two real eigenvalues and Xf) and 
two real eigenvectors can be found. The tensor in 
2D has 3 independent components. These 3 parameters 
determine 2 eigenvalues and the rotation angle that de¬ 
scribes orientation of the orthogonal eigenvectors with 
respect to the reference coordinate system. Let us asso¬ 
ciate with each atom i an ellipse with principal axes ori¬ 
ented along the eigenvectors of and having lengths 
Xj and Xf, as depicted in Fig. 1. 

Previously atomic level stresses were discussed mostly 
in 3D. In 3D symmetric atomic level stress tensors have 6 
independent components. Thus, previously, in particular 
in discussions related to the atomic level stress energies, 
it was assumed that the local atomic environment of an 
atom is described by 6 independent stress components. 
However, in view of the present considerations, it is clear 
that if the atomic level stress tensor is diagonalized then 
its 3 eigenvalues describe the geometry of local atomic 
environment, while its 3 eigenvectors describe the ori¬ 
entation of the associated ellipsoid with respect to the 
chosen coordinate system. 

In model metallic glasses in 3D atoms often have 12 
or 13 nearest neighbours [1, 2, 17]. Working with atomic 
level stresses effectively reduces the richness of all possi¬ 
ble local atomic geometries to just 3 numbers. Of course, 
that is more convenient than dealing with more numbers 



FIG. 1. Atomic level stress tensor of any atom can be diago¬ 
nalized. Obtained eigenvalues, A^ and A^, can be associated 
with the lengths of the principal ellipse’s axes. The orien¬ 
tation of the ellipse with respect to the reference coordinate 
system is given by the angle, (fi, between the longest ellipse’s 
axis and the x-axis of the reference frame. 


associated, for example, with the description based on 
Voronoi indexes [1, 2]. However, it is unclear for which 
purposes it is enough to consider only 3 numbers and for 
which purposes it may not be enough. In consideration 
of the stress correlations between two atoms in 3D there 
are 12 physically relevant parameters: 6 eigenvalues (3 on 
each atom) describe the geometries of the two ellipsoids 
and 6 parameters describe orientations of the ellipsoids 
with respect to the line from one ellipsoid to another. 
A representation in a particular coordinate frame needs 
another 3 parameters that describe the orientation of the 
line from one atom to another. 

In Ref.[15, 16] correlations between the eigenvalues 
of the same atom has been considered for 2D and 3D 
Lennard-Jones liquids. There it was argued that there 
are correlations between the stress eigenvalues of the 
same atom. 

Here we are interested in the correlations between the 
stress elements of different atoms. If we want to con¬ 
sider stress correlations between two different atoms in 
2D then we associate ellipses with both atoms and con¬ 
sider the correlations between the eigenvalues and the 
orientations of the ellipses, see Fig. 2. It is clear that in 
isotropic one component liquids all physically meaningful 
pair correlation functions should depend only on distance 

Tij. 

If the local atomic level stress tensor is known in the 
reference frame then its eigenvalues and eigenvectors can 
be found. In 2D we have: 


K 

= ( 1 / 2 ) 


( 2 ) 


{ar + af) ± ^{ar-aff+Aiarr 


tan(^/’0 = = 77(7’^ - af) ■ (3) 

Further we will assume that ipi G (—7r/2,7r/2]. Physi¬ 
cally angle (pi is defined up to an integer multiple of tt as 
rotation by angle tt does not change the ellipse. 

In the orthogonal coordinate system based on the 
eigenvectors of a particular local atomic stress tensor this 
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FIG. 2. The orientations of the ellipses with respect to the 
line connecting atoms i and j are given by the angles ijjij and 
xfjji. The orientation of the line connecting atoms i and j with 
respect to the X-axis of the reference frame is given by the 
angle 9ij. 


stress tensor is diagonal with eigenvalues and on 
the diagonal. In potentials with repulsive and attractive 
parts the values of some A can be negative. The nega¬ 
tive value of A corresponds to the case when the atomic 
environment of an atom is dilated along the eigenvector 
associated with this A. Here we assume that potentials 
that we consider are purely repulsive. Such systems held 
together at some density by periodic boundary condi¬ 
tions. In such cases, both A^ and A^ are positive and we 
order them to have A^ > A^. 


III. TRANSFORMATIONS OF STRESS 
TENSORS UNDER ROTATIONS 


In this section we provide some well known facts about 
transformations of stress tensors under rotations [23]. We 
will need these facts in our further considerations. 

Let us suppose that there are A and B coordinate 
frames in 2D and that frame B is rotated with respect 
to frame A on angle 9 in the counterclockwise direction. 
The components of the stress tensor S in frame B can be 
expressed through the components of the stress tensor in 
frame A using the rotation matrix R{9): 


Sb 


R{9)SAR^{e), 


R{e) = 


cos{9) 
— sin(0) 


sin(d) 

cos(6l) 


( 4 ) 


where R"^{9) is the transpose of R{9). In terms of com¬ 
ponents (4) leads to: 

= aY [cos(6l)]^ -I- [sin(6»)]^ -b sin(2d) , (5) 

[sin(6*)]^ + crV' [cos(0)]^ — sin(2d) , (6) 

= -(1/2) [ar - <] sin(20) + aj cos(20) . (7) 


Let us now suppose that the angle between the first 
eigenvector of atom i and the x-axis of our reference 
frame is ipi. In the frame of its eigenvectors the com¬ 
ponents of the stress tensor of atom i are af^ = Xj, 
= 0. In order to find the compo¬ 


nents of the stress tensor of atom i in our reference frame 
we should “rotate” the components of the stress tensor in 
the frame of its eigenvectors on angle —(pi using (5,6,7). 
Thus we get: 

crf"^ = A,^ cos^ipi) + Xi sin^ipi) , (8) 

sin2((^,) -b A^ cos^{ip,) , (9) 

= (1/2) (A^ - A/) sin(2(^,) . (10) 

Note also that: erf® — a^y = (Af — Af) cos(2'i/y). 


IV. CORRELATION FUNCTIONS BETWEEN 
THE ELEMENTS OF ATOMIC LEVEL STRESS 
TENSORS OF DIFFERENT ATOMS 

In this section we derive the expressions for selected 
correlation functions between the atomic level stress ele¬ 
ments in terms of eigenvalues and eigenvectors of atomic 
level stress matrices. 


A. Correlation functions in the directional frame 

It is useful to start this section from an argument which 
plays a very important role in this paper. 

Let us consider a pair of atoms i and j separated by ra¬ 
dius vector rij = rj — Vi. We associate with the direction 
of Tij a directional coordinate “r^-frame” whose :E-axis 
is along r^. The notations and erpU) will be 

used for the a/? and jS components of the stress tensors 
of atoms i and j in the r^-frame. Further we consider 
the products erp {i)ap(j) in the -directional frame 
and average such products over the pairs of atom sepa¬ 
rated by radius vector = r. It is important that this 
averaging is performed over the values of the stress ten¬ 
sor components in the representation associated with the 
rij-frame. 

For the following it is necessary to realize that for 
isotropic systems of particles the averaging, 

( 11 ) 

should not depend on the direction of r, while it can 
depend on r = jrj. This is essentially what isotropicity 
means. 


B. Transformation of correlation functions under 
rotations 

Our goal in this section is to express the correlation 
functions between the stress tensor components in an ar¬ 
bitrary frame in terms of the correlation functions in the 
rij-frame introduced in the previous subsection. 

Thus, let us express the product cr®*^(i)CT®^(j) in the 
coordinate frame which is rotated on the angle —9ij with 
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respect to rij-frame in terms of stress tensor components 
in the ry-frame. 

For this we should rotate, according to (5,6,7), the 
stress tensor components of atoms i and j in the r^- 
frame on the angle —Oij and then form the products of 
the stress tensor components in the rotated frame. From 


(7) we get: 



(12) 

= (1/4) [afj^ii) 

- <(*)] [<(4) - <(4)] [sin(2%)]^ 

+(1/4) [affiz) 

- [sin(4%)] 

+(l/4)<"(*) K^(i) - <(/)] [sin(4%))] 


[cos(2%)]^ . 


The right hand side of (12) can be expanded in terms 
containing the products of the stress tensor components 
in the ry-frame. 


Let us now average (12) over the pairs of atoms i and j 
separated by = r (this fixes the value of 9ij = 6). For 
briefness and as an example let us consider a particular 
term, cr^^(i)(T®^(j) [cos(20ij)]^, that appears on the right 
hand side of (12). In performing the averaging we get: 

«"(*)<"(j)[cos(2%)]^)..^.=. (13) 

= {<ni)<!{j))r.,=r[cOs{2ef . (14) 

In the transition from (13) to (14) the [cos(26*y)]^ was 
taken out of the averaging since the averaging is per¬ 
formed for a Hxed value of Vij = r and it also means 
that the averaging is performed for a fixed value Oij = 0. 
It follows from the previous subsection (IV A) that in 
isotropic medium {o^j {S))rij=r should not depend 

on the direction of r, but can depend on r = |r|. 

Thus, in performing the averaging of the products of 
the stress tensor components in (12), as it was done in 
(13,14), it is possible to average over all pairs of atoms 
separated by = r irrespectively of the direction of r. 
It is only necessary to ensure that the values of the stress 
tensor components on the right hand side of (12) are al¬ 
ways calculated in the directional Vij -frame corresponding 
to each pair of atoms i and j. 

It follows from the above considerations that the value 
of the correlation function {i)a^y {)))at some r 

and 9 can be expressed as a linear combination of the 
correlation functions between the atomic level stress ele¬ 
ments in the ry-frame multiplied on some functions of 9. 
Note that the dependence on 9 in (12) appears in the re¬ 
sult of rotation from the r^j-frame into the frame in which 
rij forms angle 9 with the a:-axis. Thus in an isotropic 
medium the physical essence of the atomic level stress 
correlations is contained in the correlation functions as¬ 
sociated with the -frame. In an isotropic medium these 
correlations function should depend only on distance. 


C. Expressions for the selected stress correlation 
functions in terms of eigenvalues and eigenvectors in 
the Cij-directional frame 

It follows from the two previous subsections (IV A, 
IV B) that in order to find correlation functions of the 
atomic level stress components in any coordinate frame 
it is sufficient to know correlation functions in the direc¬ 
tional rij-frame. 

It is easy to express the correlation functions in the 
rij-frame in terms of eigenvalues and eigenvectors of the 
atomic level stress matrices. Let us suppose that the first 
eigenvectors of the stress matrices of atoms i and j form 
angles ipij and ipji with the direction r^, as shown in 
Fig.2. 

From (8,9) it follows that the rotation invariant atomic 
level pressure on atom i is: 

p, ^ (1/2) = (1/2) (A) + . (15) 

Correspondingly 

{PiPj)rij=r = (1/4)((A/ -f a/) (a] -I- X^))rij=r ■ (16) 
It also follows from (8,9,10) that: 

(P.<"(j))..,=. (17) 

= (1/4)((A- -f A/) (A] - sin(2V:p))r-.,=r , 

{p^ [</ (/) - O-f/ (/)] )r-.,=r ( 18 ) 

= (1/2)((A- -f Xf) (A] - X]) C0s(2V:p))r-.,=r , 

(19) 

= (1/4)((A- - Xf) (A] - Xf) [sin(2V’ij)sin(2V’p)])r,,=r , 

(20) 

= (1/2)((A- - Xf) (A] - Xf) lsm(2f;ij) cos(2f^ji)]}r,^-=r , 

( - <"(*)] [<(J) - CrffU)] }r,,=r (21) 

= ((a! - xf) (A] - Xf) [cos(2T(ij) cos(2f;ji)]}r,^.=r , 

Note that the right hand sides of (15,16,17,18,19,20,21) 
depend on the invariant parameters of the atomic level 
stress ellipses and on their rotation invariant orientations 
with respect to the direction of . Thus, in Hnding how 
(16,17,18,19,20,21) depend on r in isotropic systems it is 
possible to average over all pairs separated by r irrespec¬ 
tive of the orientation of r. This is in agreement with the 
argument from subsection (IV A) that states that corre¬ 
lation functions between the components of atomic level 
stresses in the directional = r frame should not de¬ 
pend on the direction of r. 


D. The stress correlation function 
((j‘^^(i)cr^®'(ji))i..^=T' in the arbitrary reference frame 
expressed in terms of eigenvalues and eigenvectors. 

The correlation function {a^y{i)a^y{j))r--=r in any 
fixed reference frame depends on r, i.e., on r and 9. Us- 
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ing expressions (12,13,14,19,20,21) it is straightforward 
(although a bit tedious) to obtain the following expres¬ 
sion: 

{a-y{i)a^y{j))r.,=r (22) 

= (1/8) [J^i(r) - F 2 {r) cos(40) -|- F^{r) sin(40)], 

where 

F,{r) = {F,)r^^=r (23) 

= ((A* - A/)(A] - A|) cos(2V'ij- - 2il;ji))r^^^r , 

F2ir) = {F2U,=r (24) 

= ((A* - A/)(A] - A|) cos( 2 V'ij- -h 2lpji))r,^^r , 

Fsir) = {F^U^=r (25) 

= ((Ai ~ Aj )(Aj — \j) sin(2'(/>ij -|- 2'iljji))rij=r ■ 

Note that dependence of (22) on 9 originates from (12), 
i.e., from the rotation from the directional rij-frame into 
the coordinate frame that forms angle 9 with the direc¬ 
tion of rij = r = [rcos(0),rsin(0]. Thus the dependence 
of (22) on 9 merely reflects the rotational properties of the 
stress tensors. Also note that all physically meaningful 
information about correlations between the parameters 
of atomic level stresses is contained in functions T’i(r), 
F 2 (r), and F^ir). 

In finding Fi(r), F 2 {r), and ^/(r) in isotropic medium 
the averaging can be performed over all pairs of atoms i 
and j separated by distance r irrespectively of the direc¬ 
tion of r. 

Two alternative derivations of the formulas 
(22,23,24,25) are presented in Appendices (A, B) 

In order to understand the meaning of correlation 
function Fi let us consider the contribution from some 
atoms i and j to this function. It follows from (23) that: 

1) If one of the ellipses is a circle, for example Aj = Af, 
then the contribution from this pair of atoms is zero. 
Thus correlation function Fi(r) contains contributions 
only from those pairs of atoms in which there are finite 
shear deformations of the environments of both atoms. 

2) If ellipses of atoms i and j have the same orien¬ 
tation with respect to the line connecting them then 
cos( 2 'i/y — 2ipji) = 1 and the contribution from this pair 
of ellipses is the maximum possible contribution from 
the pairs of ellipses with the same distortions. 

3) If ellipses of atoms i and j are orthogonal to each 
other, i.e., ifij = ifji ± 7r/2 then cos(2-!/)y — 2'f’ji) = —1 
and the contribution from this pair is the minimum 
possible contribution. 

4) If tpij = ifji ± 7 r /4 then the contribution is zero. 

Note also the following. If large axes of the ellipses of 
atoms i and j are aligned then these ellipses have the 
same orientation with respect to any line, not only the 
line that connects them. Thus it is likely that rather 
simple organization of ellipses provides a maximum 
to the function Fi. It is the organization when all 
ellipses have the same shear distortions and the same 


orientations. This observation might be of interest for 
understanding the nature of viscosity. It follows from the 
Green-Kubo expression that viscosity is determined by 
decay in time of the function Fi{r), i.e., for calculations 
of viscosity it is necessary to consider stress of atom i at 
time zero and stress of atom j at time t {F 2 {r) does not 
contribute since integration over 9 in (22) leads to zero). 

In order to understand the meaning of correlation 
function F 2 {r) from ( 24 ) note the following: 

1) As in the case with Fi{r), only pairs of atoms in 
which both atoms have shear distortions contribute. 

2) The maximum contribution, for the given distortions, 
comes from the ellipses for which ipij = —'ipji, i.e., from 
those ellipses whose orientations are mirror-symmetric 
with respect to the line connecting them. 

3) If the deviation from the mirror symmetry is 7 r/ 2 , i.e., 
xfij = —'ipji ± 7r/2 then the contribution is the minimum 
possible contribution. 

4) If Tpij = —'tpji ± 7 r /4 then the contribution is zero. 

Due to a mirror symmetry we must have F^ij) = 0. 
This is because reflection with respect to the direction 
from i to j changes the signs of angles tpij and f/ji) but 
does not change the eigenvalues. In our simulations 
F 3 (r) averages to zero up to the noise level. 


E. Stress correlation function 

From (7,15,17,18), similarly to how it was done for 
{a-y{i)a-yij))r.,=r, we get: 

(26) 

= (1/4) [F 4 {r) cos(2d) -|- F 5 (r) sin(2d)] , 

where 

Fiir) = ((A/ -f A2)(A] - A|) sin{2iPji))r,^^r , (27) 

F^ir) = ((A/ + Xf)iX] - X]) cos(2^,,))r.,=r • (28) 

In finding F^lr) and F^{r) in isotropic medium the av¬ 
eraging can be performed over all pairs of atoms i and j 
separated by distance r irrespectively of the direction of 
r. 

Due to mirror symmetry, the function F 4 should 
average to zero (it does in simulations). 

In order to understand the meaning of F 5 from (28) 
note the following: 

1) The larger is the pressure on atom i and the shear 
distortion of atom j, the larger is the contribution from 
this pair to F 5 . 

2) If the ellipse of atom j is aligned with the direction 
from i to j then cos(2rpji) = 1 and there is the maximum 
possible contribution for the given ellipses’ shapes. 

3) If the ellipse of atom j is orthogonal to the direction 
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from i to j then cos{2'il)ji) = — 1 and there is the minimum 
possible contribution for the given ellipses’ shapes. 

4) If ipji = 7r/4 then the contribution is zero. 


F. Stress correlation function (((rf’” — cr^^)a^^)rij=r 

From (5,6,7) and (19,20,21), similarly to how it was 
done for we get: 

( 29 ) 

= (l/2)Fi(r) sin(40) + Fe{r) cos{4:9) , 

where Fi(r) is given by expression (23) and: 

Feir) = ((A- - A,^)(A] - A|) cos(2V’*j) sin(27/>ji))r,^-=.r.(30) 

In finding Fq (r) in an isotropic medium the averaging can 
be performed over all pairs of atoms i and j separated 
by distance r irrespective of the direction of r. 

The function Fq should average to zero due to mir¬ 
ror symmetry with respect to the direction from i to j 
since under reflection cos{'ipij) does not change sign, while 
sin(ipji) does. We verified this in our simulations. 


G. Simpler correlation functions and normalization 
of the correlation functions 

Correlation functions 7^1,2, 3 , 4 , 5 ,6 are somewhat compli¬ 
cated as they represent averages over three or four param¬ 
eters. Before considering them it makes sense to consider 
simpler correlation functions which represent averaged 
products on two parameters only. It is expectable that 
stresses of particles which are far away from each other 
are not correlated. This makes it reasonable to consider 
the following correlation functions: 

G,,{r) = {l/Zlml + A2)(A] + A])}r,,=r - 1 , (31) 

G^™(r) = (l/z!)((Ai-A2)(A]-A2)),.^.=,.-l , (32) 
G™p(r) = (l/Z+Z.){(Xl - A2)(A] + X])}r„=r - 1 (33) 
C2±ir) = {cos{2'ipij ± 2ll)ji))ri^=r , (34) 

where Z± = (A) ± Af). 

Functions Gpp{r), Gmm(,r), and Gmpir) describe cor¬ 
relations between the eigenvalues (or eigenstresses) of the 
stress matrices of atoms i and j without taking into ac¬ 
count the orientations of the eigenvectors. Note that 
since Pi = (1/2)(A^-|-A^) the function Gpp{r) from (31) is 
directly related to the pressure-pressure correlation func¬ 
tion between atoms i and j. It follows from Appendix A 
and formula (2) that the function Gmm represents cor¬ 
relations between the total amounts of shear on atoms 
i and j. Finally, Gmp{r) describes correlations between 
the total shear on atom i and the total pressure on atom 
j. Functions G 2 ±(r) from (34) describe correlations in 
the orientations of the eigenvectors of the stress matrices 


of atoms i and j without taking into account the magni¬ 
tudes of the eigenvalues. 

It is also reasonable to introduce normalized versions 
of the correlation functions Ai 2 , 3 , 4 , 5 , 6 : 

Fl,2,Zfi{Tij) = Ai,2,3.6/^- ) (35) 

A4,5 = F^^^I{Z+Z_) . (36) 

V. ANALOGY WITH THE ESHELBY’S 
INCLUSION PROBLEM 

In this section we discuss from the perspective of the 
Eshelby’s inclusion problem [21-25] the stress correlation 
function which is analogous to the atomic level stress cor¬ 
relation function discussed in the previous 

section. In drawing this analogy it is assumed that the 
central atom i is analogous to the Eshelby’s inclusion (J) 
that generates a stress field in the matrix at point J (on 
atom j). In particular, we argue that the angular de¬ 
pendence of the {(T^^(T^^)r=r stress correlation function 
obtained in Ref. [14] is related to the rotational proper¬ 
ties of the stress tensors and not to the anisotropy of the 
stress Held associated with the Eshelby’s solution. 

There are two points which we need from the the Es¬ 
helby’s solution. 1) The final strain and stress fields in 
the inclusion after the deformation, placing the inclusion 
back into the matrix, and joining are constant. The final 
strain and stress fields in the inclusion, of course, depend 
on the unconstrained strain initially applied to the inclu¬ 
sion. 2) If we know the unconstrained strain applied to 
the inclusion then the final strain and stress fields in the 
inclusion and in the matrix can be found. Further we as¬ 
sume that there is a one-to-one correspondence between 
the stress fields in the inclusion and in the matrix. See 
also Appendix C. 

We are interested in the correlation functions between 
the inclusion (/) and some point (J) in the matrix. Sim¬ 
ilarly to how it was done for the atomic level stresses, we 
can associate with I the stress ellipse whose parameters, 
(A}, Aj), and whose orientation, ipij, with respect to rjj 
are known. The fact that the stress field is the same 
everywhere in the inclusion serves well for this purpose. 

Since the inclusion’s stress ellipse is known, the stress 
field at any point J in the matrix can be found. Since 
the stress tensor at point J is known it can be diago¬ 
nalized and thus it is possible to associate with point J 
its own stress ellipse with parameters (Aj, Aj) and the 
orientation 'ipjj with respect to r/j. 

At this point it becomes apparent that considerations 
of correlations for the Eshelby’s inclusion problem are 
quite similar to the considerations that were already done 
for the atomic level stresses. There is, however, an impor¬ 
tant difference. Thus, in the case of atomic level stresses 
correlations between the parameters of atomic level stress 
ellipses have a probabilistic character. In contrast, in the 
case of the Eshelby’s inclusion problem the stress field in 
the inclusion deterministically defines the stress field in 



7 


the matrix. Thus, Aj, Aj, and ipJi are the functions of 
A}, Xj, i^ij, and r/j: 

(^/) ^ 7 ) V'/J) ) (37) 

(^ 7 ; ^ 7 ) V'/Jj ^7j) ) (38) 

^J7 = V’J7 (a}, Aj,V'7J,?"7j) • (39) 

In (37,38,39) angles ipu and ipji are the angles between 
the larger ellipses’ axes and the direction rjj. Note that 
in isotropic elastic medium Aj, Aj, and ipji should not 
depend on the direction of rjj (the direction of the in¬ 
clusion’s deformation with respect to rjj is taken into 
account by the angle ipij). 

Note that the properties of the Eshelby’s solution are 
embedded into (37,38,39). These functions, in our view, 
represent the essence of the Eshelby’s solution. In Ap¬ 
pendix C a particular case of the inclusion’s shear trans¬ 
formation is discussed and functions (37,38,39) for this 
case are derived. 

Expressions for the stress correlation functions be¬ 
tween the inclusion and the matrix can be derived in the 
same way as the expressions (22,23,24,25,26,29) for the 
atomic level stress correlation functions. For the product 
for example, we get: 

a^y(I)a^y{J) (40) 

= (1/8) [T! - J-| cos(407j) + J"! sin(4d7j)], 

where 

jre = (^1 _ _ ^2 ) cOs(2V:7J " , (41) 

.FI = (A/ - Xj){X) - X^j) cos(27/7j + 27/j7) , (42) 

•^3 = (^7 ~ ^ 7 )(^j ~ ^j) sin(27/>7j -I- 2ipji) . (43) 

The upper index e in the formulas above originates from 

the word “elastic”. Note again that Aj, Aj, and ipji in 
(41,42,43) are the functions of A}, Aj, ipu, and rjj. Also 
note that F"/, F"!, F| do not depend on 9ij. Thus in 
(41,42,43) 

= T^{X],Xj,'ipi,riA) , where n = 1,2,3, (44) 

i.e., functions Ff, F|, and F| are determined by how the 
stress field in the inclusion determines the stress field at 
point J. 

Now we comment on the connection between the func¬ 
tions Ff, F|, and F| from (41,42,43) and the functions 
Fi, F 2 , and F 3 from (23,24,25). The functions Ff, F|, 
and F| are written for a particular set of values A}, Xj, 
ipij, and rjj. In order to draw a parallel with the atomic 
level stress correlation functions in liquids it is necessary 
to average the functions F|, F|, and F| over the possible 
values of A}, Aj, and ipu which can be associated with 
the parameters of the inclusion’s stress ellipse. Thus: 

Kirij) = (F|(A},A?,V'7,r7j))AbA?,^.,, , (45) 

where n = 1,2,3. In (45) it is presumed that every set 
of parameters at I deterministically leads to certain pa¬ 


rameters at J via the Eshelby’s solution. In (45) there is 
no averaging over the distance (scalar) rjj. Correspond¬ 
ingly functions F| depend only on rjj = r. 

In liquids there is no deterministic relation between 
the parameters and orientations of the atomic level stress 
ellipses of atoms i and j. In liquids there is only a proba¬ 
bilistic relation. Thus in calculations of Fi(r), F 2 (r), and 
F;^{r) in liquids (23,24,25) the averaging goes not only 
over A), Xf, ipij, but also over Aj, Aj, ipji. Implicitly in 
calculations of (23,24,25) there is also the averaging over 
the directions of for a fixed value of r^ . Since it is 
assumed that the undistorted inclusion and the matrix 
are isotropic there is no need to average (45) over the 
directions of rjj. 

Note that if {cr ^were calculated from (B2) in a 
particular reference frame, by averaging over the possi¬ 
ble distortions of the inclusion, it still would depend on 
9ij. This dependence, however, would not reflect the 
essence of the angular dependent Eshelby’s stress field. 
The dependence on 9jj in (B2) reflects the rotational 
properties of the stress tensor. The angular dependen¬ 
cies observed in Ref. [14] correspond to the dependence of 
{a“^y{I)a“^y{J)) on 9ij in (B2). This is not the angular 
dependence of the Eshelby’s field. The angular depen¬ 
dence of the Eshelby’s stress field is embedded in how 
Aj, Aj, iiji depend on A}, Aj, and rij. 

VI. RESULTS OF MD SIMULATION 
A. Stress correlation functions 

In our Molecular Dynamics (MD) simulations we con¬ 
sidered the same 2D system of particles that has been 
studied in Ref. [14]. We used the same Yukawa poten¬ 
tial and the same LAMMPS MD program [18, 19]. We 
studied the systems of two sizes. In the small system 
the number of particles was N = 2500, while the dimen¬ 
sions of the rectangular periodic box were = 50.1021, 
Ly = 43.3897. Our small system has the same size as 
the system studied in Ref. [14]. Another (large) system 
that we studied contained N = 22500 particles, i.e., nine 
times more than the small system. The dimensions of 
the large system were = 150.306, Ly = 130.169. The 
particles’ number densities in the small and large sys¬ 
tems are the same. We performed simulations in NVT 
and NVE ensembles. 

In all cases the systems were prepared by melting tri¬ 
angular lattice at reduced temperature T = 5. After the 
equilibration at T = 5 the temperature was reduced in 
several steps that followed by equilibration at every tem¬ 
perature (in NVT ensemble) or at every value of fixed 
total energy (in NVE ensemble). The temperature in 
NVT ensemble was introduced via Nose-Hoover thermo¬ 
stat. The damping parameter corresponded to 100 MD 
steps and also to 0.1 of the time unit. 

In our simulations, we reproduced the dependence of 
potential energy on temperature presented in Fig. 1 of 






FIG. 3. In every panel the curves from the top to the bottom correspond to temperatures: T = 3, T = 2, T = 1.4, and T = 1. 
The curves were shifted vertically for the clarity of the presentation, (a) Pair density function, (b) Gmm correlation function 
from (32). (c) Gmp correlation function from (33). (d) Gpp correlation function from (31). 


Ref. [14]. 

Atomic configurations for calculations of the correla¬ 
tion functions related to the eigenvalues and eigenvec¬ 
tors of atomic level stresses were collected on the small 
system in the NVE ensemble at total energies which 
corresponded to the following temperatures: T(3) = 
3.06 ±0.04, T(2) = 1.97 ±0.03, T(1.4) = 1.43 ± 0.02, 
r(l) = 0.99 ± 0.02. The averaging was done over 1000 
configurations at every temperature. For the tempera¬ 
ture T = 1 the time interval between the two consecutive 
configurations was 10^ MD steps. Each MD step cor¬ 
responded to 0.001 of the time unit. During these 10^ 
MD steps the mean square atomic displacement reaches 
~ 1.38(7. 

Different correlation functions per pair of particles are 
shown in Figs.3,4,5. The dependencies of the functions 
Fi, ±2 and F 5 , i.e., all the non-zero ones, on distance are 
shown in Fig. 5. At T = 1 we have Z_ = (A^ — A^) « 
10.82 and Z+ = (A^ ± Xf) « 40.48. 

Figures 3,4,5 demonstrate that there are r^-dependent 
correlations between the parameters of the atomic level 
stress ellipses and in their orientations. These correla¬ 
tions gradually decrease with increase of r^-. It is clear 
that functions Gmmir) in Fig.3(b) and < cos(2(/>jj — 


2V’ji) > in Fig.4(a) exhibit more pronounced changes 
than does PDF [Fig.3(a)] on decrease of temperature. It 
is also clear that the first peaks in Fi and ±2 [Fig.5(a,b)] 
also demonstrate more pronounced changes on decrease 
of temperature than does PDF. However, it is also more 
difficult to interpret these changes. Yet, developing fea¬ 
tures in {cos{2tpij — 2 i/)ji)) suggest that some ordering 
happens in the mutual orientations of the ellipses as¬ 
sociated with the atoms separated by the distance cor¬ 
responding the first minimum in the PDF. There also 
appears to be a certain similarity in the behaviours of 
(cos(2(/:ij — 2'ijjji)) and Fi. This similarity suggests that 
changes in Fi are caused by changes in {cos{2'tfjij — 2tpji)). 
See expression (23) for Fi. Thus changes in Fi are likely 
to be caused not by changes in the eigenvalues of the 
stress ellipses, but by changes in the mutual orientations 
of the ellipses. However, also note that there are changes 
in Gmmir) in Fig.3(b). 

Figure ( 6 ) shows how the function Fi{r)p{r) changes 
with temperature. It follows from the figure that as tem¬ 
perature is reduced there develops a pronounced min¬ 
imum at the position of the first minimum, rG^, of 
p{r). Thus changes in Fi{r) are also well observable 
in Fi{r)p{r) despite the fact that the number of atomic 
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FIG. 4. Evolutions with temperature of the functions 
{cos{2'tpij — 2tpji)) and {cos{2^ij + 2tpji)). In both panels 
the curves from the top to the bottom correspond to temper¬ 
atures T = 3, T = 2, T = 1.4, and T = 1. The curves were 
shifted vertically for the clarity of the presentation. 


pairs separated by is relatively small. 

The curves in Fig.5 can be converted into the 2D in¬ 
tensity plots equivalent to those presented in Ref. [14] us¬ 
ing formulas (22,26,29). Thus, if we want to find the 
stress correlation function at a point with coordinates 
(x, y) we define r = y/x"^ + y'^ and 9 = arctan (y/x). Us¬ 
ing these values in (22,26,29) the stress field at {x,y) 
can be found. This conversion applies because for 
particles i and j with coordinates {xi,yi) and {xj^yj) 
the values of Vij and 9ij that go into the formulas 
(22,26,29) are = y/{xj - Xi)'^ + {yj - yi)'^ and % = 
arctan((j/j — yi)/(xj — Xi)). However, in making the 21) 
stress correlation function plots it is assumed that the 
particle i is at the origin. 

The results of the conversion described above for T = 1 
are presented in Fig.7. It is obvious that the 2D plots in 
Fig.7. are very similar to those shown in Fig.5 of Ref. [14]. 
Note that the 2D plots presented in Fig.7 were obtained 
from only 3 functions, i.e., Fi(r), F 2 (r), and ^^(r) which 
depend only on r. This proves that the dependencies on 
9 presented in the 2D plots in Ref. [14] follow from the 
tensorial rotational properties. 



FIG. 5. Evolutions with temperature of the normalized func¬ 
tions Fi from (23,35), F 2 from (24,35), and F5 from (28,36). 
In every panel the curves from the top to the bottom corre¬ 
spond to temperatures T = 3, T = 2, T = 1.4, and T = 1. 
The curves were shifted vertically for the clarity of the pre¬ 
sentation. 



FIG. 6. Evolution with temperature of the function Fi (r)p(r). 
The curves from the top to the bottom correspond to tem¬ 
peratures T — 3, T = 2, T = 1.4, and T = 1. The curves 
were shifted vertically for clarity of presentation. The dashed 
curve shows the scaled p{r) at T — 1.0. As temperature 
decreases the first minimum in p(r) becomes deeper. This 
deepening overlaps with the development of the minimum in 
Fi{r). Thus, the development of the feature in Fi{r) at the 
position of the first minimum of p(r) is also well pronounced 
in Fi{r)p{r) . 
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FIG. 7. (a) Atomic level stress correlation function See formula (22). The functions Fi and F 2 in were 

normalized according to (35). The function F 3 is zero, besides the noise. It is clear that panel (a) of this figure is very similar 
to the panel (a) of Fig.5 in Ref.[14]. (b) Atomic level stress correlation function (piuj^). See formula (26). The function F 4 in 
{piO-j^} averages to zero. Thus only function F 5 is left. In order to produce this hgure we had to subtract the average pressure 
from the diagonal components of the atomic level stress tensor. Effectively this means that the value of A averaged over A^ 
and A^ of all atoms, i.e., A^"®, was subtracted from the values of A^ and A^ of every atom. Because of this subtraction we can 
not use normalization (33) since (A^ + A^ — 2A“”®) averages to zero. Thus we used normalization (32) instead. We also scaled 
intensity on the z-axis by a factor of 4. It is obvious that panel (b) of this figure is very similar to the panel (b) of Fig.5 in 
Ref.[14]. (c) Atomic level stress correlation function — o'f“)crJ^). See Eq.(29). Only function Fi in {{af^ — cr™)(Tj^) is 

non-zero. We scaled the function by a factor of 10 along the 2 -intensity axis. It is clear that correlation function in panel (c) 
is rather similar to the correlation function in panel (c) of Fig.5 in Ref.[14]. There is a difference with the two “circles” at 
{f/O'} ~ 1- It is possible that authors of Ref. [14] did not look into such small distances and thus in their figure these “circles” 
fall into the central green region. Besides this difference their figures and ours look rather similar. 


B. Is system studied a true liquid or is it in a 
hexatic phase? 

Finally, we comment on the following statement made 
in Ref. [14]. It is stated there that at T = 1 the system is 
in a true liquid state, while at T = 0.95 the system is in 
a hexatic state. 

In order to make a distinction between the true liq¬ 
uid and haxatic states it was assumed in Ref. [14] that 
in a true liquid state bond-order correlation function de¬ 
cays exponentially with increase of distance, while in the 
haxatic state the bond-order correlation function decays 
algebraically. 

We calculated how the bond-order correlation function 
depends on distance in systems of two sizes. In the small 
system containing N = 2500 particles {Lyj2) = 21.7, 
while {Lyj2) = 65.1 in the large system with N = 22500. 
The small system was used in Ref. [14]. In our cal¬ 
culations we assumed that two atoms are the nearest 
neighbours if they are separated by a distance smaller 
than the position of the first minimum in the PDF, i.e., 
{rijja) < 1.36cr. The results are presented in Fig.8,9. 

It follows from Fig.8,9 that on decrease of temperature 
haxatic order undoubtedly develops in the systems. The 
comparison of Fig.8 with Fig.9 suggests that at T = 0.95 
the small system exhibits observable size effects. Note 
that at T = 0.95 the decay length is larger than (1/2) 
of Lyl2 in the small system. Thus in the small system 
the bond-order correlation function does not decay com¬ 


pletely on the length of the half of the simulation box. It 
also follows from the data obtained on the large system 
that exponential fit to the data is better than can be any 
algebraic fit at both temperatures. Thus, in our view, 
it follows from Fig. 8,9 that it is impossible to make a 
qualitative distinction between the liquid states at T = 1 
and T = 0.95. The observation of the algebraic decay at 
T = 0.95 reported in Ref. [14] is probably related to the 
size effects. 


VII. CONCLUSION 

It was demonstrated that it is possible to study liquid 
(and glass) structures through considerations of corre¬ 
lations between the eigenvalues and eigenvectors of the 
atomic level stress tensors of different atoms. It was 
shown that on decrease of temperature some of the stud¬ 
ied correlation functions exhibit pronounced changes in 
the range of distances that corresponds to the first min¬ 
imum of the pair density function. These changes could 
not be guessed from the behaviour of the pair density 
function. Thus the suggested method provides additional 
information and it is of interest to investigate evolution 
of stress correlations with this method in model super¬ 
cooled liquids on decrease of temperature. 

We also demonstrated that interpretations of the an¬ 
gular dependencies of the stress correlation functions re¬ 
ported in Ref. [14] are essentially incorrect. In particular. 
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FIG. 8. Bond-order correlation functions for the system of 
N = 2500 particles. There are 10 cyan curves in the figure. 
Every cyan curve is the average over 1000 independent config¬ 
urations. The mean square displacement between the consec¬ 
utive configurations was larger than the interatomic distance. 
There are three blue curves in the figure. They represent the 
mean over the 10 cyan curves and the average ± the error 
of the mean. There are 10 green curves in the figure. Every 
green curve represents the average over 100 independent con¬ 
figurations. There are 3 magenta curves which represent the 
average over the green curves and the average ± the error of 
the mean. The red and maroon curves show the fits obtained 
from the larger system of = 22500. 



FIG. 9. Bond-order correlation functions for the system of 
N = 22500 particles. There are 21 cyan curves in the figure. 
Every cyan curve is the average over 100 independent config¬ 
urations. The mean square displacement between the consec¬ 
utive configurations was larger than the interatomic distance. 
There are three blue curves in the figure. These blue curves 
represent the mean over the 21 cyan curves and the mean ± 
the error of the mean curves. It is possible that at distances 
larger than r/cr = 30 at T = 0.95 there again appear size ef¬ 
fects. There are three green curves in the figure. Every green 
curve is the average over 200 independent configurations. The 
red and maroon curves are the fits to the data. 


the authors of Ref. [14] associate the angular dependen¬ 
cies observed in the stress correlation functions with the 
angular dependencies of the Eshelby’s stress field. We 
demonstrated that anisotropic stress fields observed in 
Ref. [14] originate from the rotational properties of the 
stress tensors. We also had shown that information which 
is really related to the anisotropic Eshelby’s stress fields 
is embedded into the isotropic stress correlation functions 
E’i(r), F 2 {r), and E 5 (r) which we studied in this work. 

From a purely pragmatic perspective we have shown 
that eight 2D-panels of the stress correlation functions 
presented in Ref. [14] can be reproduced using only 3 cor¬ 
relation functions which depend only on r, i.e., from the 
function Fi{r), F 2 {r), and E 5 (r) [20]. This clearly ad¬ 
vances understanding of the plots of the stress correlation 
functions presented in Ref. [14]. It also follows from our 
results that instead of studying distance dependence of 
the integrals of the 2D stress correlation functions over 
some angles, as it has been done in Ref. [14], it is more rea¬ 
sonable to study how functions Fi{r), F 2 {r), and E 5 (r) 
depend on distance. 

We also demonstrated that because of size effects the 
distinction made in Ref. [14] between the normal liquid 
and haxatic states is invalid. 


Appendix A: Alternative derivation of 
structure 


In 2D in a particular reference coordinate frame nu¬ 
merical representation of the atomic level stress tensor 
(7 is a 2 X 2 matrix. This matrix is real and symmetric 
{i.e., = (y^^), thus it can be diagonalized. We can 

work directly with its components or with corre¬ 

sponding pressure p and two shear components, Si and 
S 2 ; or with real eigenvalues Ai ^2 and the orientation of 
two orthogonal eigenvectors: 




P+Si S2 

a^y ayy 


S2 p-Si_ 


= R{<p) 


Ai 0 
0 A 2 


{R{p)f. 


Here R{ip) is the 2x2 matrix of rotation in positive (or 
counterclockwise) direction by angle p: 


R{p) 


cos (p — sin if 
sin If cos tp 


Pressure and shear components are expressed through 
cr“^ components as 

p= l{a^^+ ayy) , Si = S 2 = . 

It will be convenient for us to combine the shear com- 
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ponents into a single complex number s = si + z S 2 - The 
total amount of shear is given by its absolute value: 

kl = \lsi + si = ^\{a^^ - avyy + , 

while the argument of s is related to the shear’s direction. 

Consider three reference frames {x,y), (x',y'), and 
{x",y"), with {x',y') being obtained from {x,y) by a ro¬ 
tation in negative (clockwise) direction by angle a, while 
{x,y) and {x",y") are mirror reflections of each other 
with respect to x-axis. We will write down the quanti¬ 
ties in (x', y') and [x", y") frames with prime and double 
prime symbols, respectively. The transformation proper¬ 
ties of the stress tensor (5,6,7) result in 

s' = s exp(2z a ), s" = s* , 

where •* denotes complex conjugation. 

Since a^y = (s — s*)/2z, we have: 

(afVf) = i [Fi(r, 9) + F*(r, 9) - F^r, 9) - F*{r, 0)] , 

(Al) 

where 

Fiir,9) = 4(s,s*), F2{r,9) = 4(s,s,); (A2) 

All the averages (•) are taken over the pairs of atoms i 
and j with = r and 9y = 9. 

By checking how s and the angle 9 are transformed by 
rotations (s' = s exp(2z a) and 9' = 9 + a) vie get 

Fi{r,9 + a) = Fi{r,9), 

F 2 (r, 0-I-a) = F 2 (r, 0) exp(4z a). (A3) 

The function Fi{r,9) does not depend on the angle 9 
at all. By considering F^ (r, 9) we exchange the roles of 
atoms z and j, which is equivalent to the change 9 —)■ 
9 + 71. Thus Fl(r,9) = Fi{r,9 + tt) = Fi{r,9), i.e., we 
get Fi = Fy All this means that Fi (r, 9) = Fi (r) is a 
real function of a single parameter r. 

If we put 0 = 0 in (A3), we get F 2 (r, a) = 
F2{r,0) exp{4:ia). Mirror reflection (s" = s* and 
9" = —9) leads to F2{r,—9) = fKr,0). In partic¬ 
ular, F 2 {r) = F 2 {r, 0 ) = F 2 {r, 0 ) is also a real func¬ 
tion of just the distance between the atoms r. Also, 
F2{r,9) = F 2 {r) exp(4z0). 

Putting these results for Fi and F 2 into the expression 
(AI) we finally get 

(crf^cr^) = I (Ai(r) - F 2 {r) cos(46l)) . (A4) 

Note that the angular dependence of this correlation 
function was obtained solely by checking how the atomic 
level stress tensors are transformed under rotations (and 
mirror reflections). Thus the physical properties of the 
liquid prescribe the r-dependence of Fi 2 {r), but not the 
0-dependence in (A4). 


Appendix B: Yet, another derivation of the 
expression for 

Let us suppose that the first eigenvectors of atoms z 
and j form angles (pi and tpj with the x-axis of our ref¬ 
erence coordinate frame. See Fig.2. 

From (10) the product in our reference coordi¬ 

nate frame has the form: 

= (l/4)(Al - - A2)sin(2(^,)sin(2(^,) . (Bl) 

Note that the dependence on angles (pi and (pj appears 
in (Bl) from the “rotations” (10) of the stresses from the 
coordinate frames of their eigenvectors into our reference 
coordinate frame. Thus dependence of (Bl) on cpi and (pj 
reflects transformational properties of the stress tensors 
under rotations. 

We then, using Fig.2, express angles (pi and p>j through 
the angles tpij, ipji, and 0^-: Lpi = ipij + 9ij and Lpj = 
ifiji -I- 9ij. Substitution of these expressions for Lpi and ipj 
into (Bl) (with the following averaging) leads to: 

(crf*'crj*')e.^.=e = (1/8) [Fi - F 2 cos(40) -f F 3 sin(40)] ,(B2) 

where Fi, F 2 , and F 3 are given by expressions (23,24,25). 
Note that (B2) is identical to (22) and (A4). 

Since the dependence of (Bl) on (pi and Lpj appeared 
from the rotational properties of the stress tensors the 
dependence of (B2) on 0 also reflects the rotational prop¬ 
erties of the stress tensors. 


Appendix C: Eshelby’s stress field in the directional 
frame for a case of shear deformation of a circular 
inclusion. Functions FI, T 2 , 

In this section we derive the expressions relating the 
eigenvalues and eigenvectors of the stress fields in the in¬ 
clusion and in the matrix for a particular case of uncon¬ 
strained shear strain applied to the inclusion. Then we 
calculate functions F^, JTf, F^ for the considered exam¬ 
ple. We start from the known formulas for the Eshelby’s 
stress field [21-25]. In particular, we use the expressions 
provided in Ref. [25]. 

We consider a particular case of unconstrained shear 
strain applied to the initially circular inclusion: 

4/3 = e* (2han/3 - 4/3) , (Cl) 

where n is a 2-dimensional unit vector that determines 
the “direction” of deformation: 

4 = cos(z/> 7 j) , fly = sin(z//j) . (C2) 

The expression for the final stress field in the inclusion 
(in the absence of external driving force) from formula 
(14) of Ref. [25] is: 

^a/3 ~ 9^aP ; 9 = ’ (C3) 
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where E is the Young’s modulus, while v is the Poisson’s 
ratio. The eigenvalues and eigenvectors of the stress ten¬ 
sor (C3) can be easily found: 

A} = +ge*, V} = [ -hcos('!/'/j), -hsin(V)/j) ], (C4) 
Aj = -ge*, Vf = [- sin(V’7j), + cos(V'/j) ] • (C5) 

The expression for the final stress field in the matrix, 
according to formula (A25) of Ref. [25], is: 

where 

(C7) 

=-A{l/rf{{l-2v) + {llrf} 

■ {(nr) {flafp + flpfa) - fafp} 

+(l/f)2 {2(1 - 2u) + [l/ff] {2n^h0 - 

-4(l/f)^ (1 - 2(1/ff] |2 (nrf - l| f^rp 

+4(l/rT{l-(l/rT} 

• |(nf) (haTp + hpfa) -2(nrf Tafp^ 

+2(l/rf (1 - (l/ff) |2 (nrf - l| . 

We are interested in the expression for the stress field in 
the coordinate frame associated with the direction from 
/ to J. In this frame r = (1,0), while (nr) = cos(ipij). 
Also note that fixfx = cos(i/'/j) and nyTx = sin('ipjj), 
while hxfy = 0 and fiyfy = 0. It is straightforward to 
obtain from (C7) the following expressions: 

{2-3(“) |sin(27/-7j), (C8) 

[•■4a. = (^) |-4(l-4+3(“) I cos(27/:/j) , (C9) 

= ( 7 ) }cos(2V’/j). (CIO) 

Using expressions (C6) and (C8,C9,C10) for the stress 
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